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Abstract 

The Barut-GirardeUo coherent states (BG CS) representation is extended to the non- 
compact algebras u{p, q) and sp{N, R) in (reducible) quadratic boson realizations. 
The sp{N, R) BG CS take the form of multimode ordinary Schrodinger cat states. 
Macroscopic superpositions of 2"~^ sp{N, R) CS (2" canonical CS, n = 1, 2, . . .) are 
pointed out which are overcomplete in the iV-mode Hilbert space and the relation 
between the canonical CS and the u{p, q) BG-type CS representations is established. 
The sets of u(p, q) and sp{N, R) BG CS and their discrete superpositions contain 
many states studied in quantum optics (even and odd A^-mode CS, pair CS) and 
provide an approach to quadrature squeezing, alternative to that of intelligent states. 
New subsets of weakly and strongly nonclassical states are pointed out and their sta- 
tistical properties (first- and second-order squeezing, photon number distributions) 
are discussed. For specific values of the angle parameters and small amplitude of 
the canonical CS components these states approaches multimode Fock states with 
one, two or three bosons/photons. It is shown that eigenstates of a squared non- 
Hermitian operator A"^ (generalized cat states) can exhibit squeezing of the quadra- 
tures of A. 



1 Introduction 

In the recent years there has been much interest in applications and generalizations of 
the Barut-Girardello (BG) coherent states (CS) § g |, g |, |, 0]. The BG CS were 
introduced ^ as eigenstates of the lowering Weyl operator K_ of the algebra sn(l, 1). The 
BG CS representation has been used for explicit construction of squeezed states (SS) for 
the generators of the group SU{1, 1) which minimize the Schrodinger uncertainty relation 
for two observables |Q and of eigenstates of general element of the complexified algebra 
su'-^ {1,1) 0, The overcomplete families of eigenstates of elements of a Lie algebra 
were called algebraic CS and algebra eigenstates 0, ||. The idea to construct SS for 
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quadratures of any non-Hermitian operator A as eigenstates of complex combinations 
uA + vA^ was put forward in ref. 0, where such eigenstates |z, m, v) were constructed for 
A = J_ and A = K_, J_ and K_ being the Weyl lowering operators of su{2) and su{l, 1) 
Lie algebras correspondingly. The su{l, 1) BG CS differ from the SU{1, 1) group related 
CS (see |T^ and references therein): the highest weight vector is the only one common 



state, while the su{l, 1) algebra related CS \z,u,v; k) of ref. [|l| contain the whole set of 
SU{1, 1) group related CS with symmetry. The general set of algebra related CS always 
contains the corresponding group related CS with symmetry as a subset. 

Passing to other algebras it is initially important to construct the eigenstates of Weyl 
lowering operators, which is a direct extension of the BG definition of sm(1, 1) CS to 
the desired algebra. The aim of the present work is to construct BG type CS for the 
symplectic algebra sp{N,R) and its subalgebras u{p,q), p + q = N, in the quadratic 
boson representation. This sp{N, R) representation is of importance in various field of 
physics |Tl|, 0. Here N is the dimension of Cartan subalgebra, while the dimension 



of sp{N, R) is A^(2A^ + l)/2, N = 1,2,...,^. 



We establish that the sp{N, R) BG CS in quadratic boson representation take the 
form of superpositions of two multimode canonical CS |jTO| \d) and \ —d) (equation (|20|) ). 
A subset of these states is found which is overcomplete in the whole Hilbert space Ti of 
the iV-mode system. Recall that the corresponding Sp{N, R) group related CS are not 
overcomplete in Ti since the representation is reducible. This property is a particular 
case of a quite general result of the overcompleteness of eigenstates of powers A^" of non- 
Hermitian Aj, j = 1, 2, . . ., provided the eigenstates \z), z = {zi, . . . , zn), of all Aj are 
overcomplete with respect to a measure independent of the phases of Zj (section 3 and 
appendix A. 2). 

Macroscopic superpositions of two canonical CS are called (ordinary) Schrodinger cat 
states Ijl^, |15]. The set of the sp{N,R) BG CS includes several subsets of ordinary 
cat states, which are extensively studied in quantum optics (see [0, |TB| and references 



therein). We introduce multimode squared amplitude Schrodinger cat states as macro- 
scopic superpositions of two sp{N, R) BG CS. Unlike the ordinary cat states these su- 
perpositions, which eventually become combinations of four mode canonical CS, can 
exhibit amplitude and squared amplitude quadrature squeezing (first- and second-order 
quadrature squeezing or hnear and quadratic squeezing) |T^, and other nonclassical prop- 
erties. Families of weakly and strongly nonclassical [l^ cat states are pointed out as 
macroscopic superpositions of two Sp{N,R) CS. There are states in these famihes that 
tend to multimode Fock states with 0, 1, 2 or 3 photons as the amplitude of their canonical 
CS components approaches zero. We note that, unlike the case of Robertson (Schrodinger) 
intelligent states [|l], |[, the cat state squeezing can not be arbitrarily strong. 

Recently [0] the BG type CS have been constructed for the u{N—l, 1) algebra. Here we 
construct overcomplete families of states for u{p, q), p + q = N, and the related resolution 
unity measures as well. We show that the 'pair CS' 0] are in fact the u(l, 1) BG CS 
\z] k) for k = 1/2, 1 . . ., while the 'two-mode Schrodinger cat states' of are particular 
case of our u{p,q) multimode squared amplitude cat states (0). 

The paper is organized as follows. In section 2 a concise review of the properties of BG 
CS representation and its relations to the canonical one- and two-mode CS representation 
(or Fock-Bargman representation) |TD| is given. An explicit relation between the two- 
mode canonical CS and the BG CS representations is obtained. Using this relation on 
can more easily establish the coincidence between the generalized intelligent states (IS) 



2 



\z,u,v; k) m and many other one and two-mode states, constructed by other authors as 
eigenstates of ua"^ + fa^^ or uab + fa^6^ ^ For example, for real u, v the states 



\z, u, v; k = l/2) coincide with the 'pair excitation-deexcitation CS' |£T|, while for real u, v 
and /c = (1 + \q\)/2 they are identical to the 'two-mode intelligent SU{1, 1) CS' p2[| . 

In section 3 the BG CS are extended explicitly to the algebra sp{N, R) in the (re- 
ducible) quadratic boson representation and overcomplete in whole H families of such 
states are constructed. Overcomplete families of eigenstates of the power 2n of Weyl op- 
erators ajOj are also built up. These states take the form of macroscopic superpositions 
with 2" canonical CS components. The u{p, q) BG type CS are considered in section 4 
(overcompleteness, resolution unity measure, particular cases and relation of their ana- 
lytic representation to that of canonical CS). In section 5 the statistical properties (weak 
and strong nonclassicality, amplitude and squared amplitude quadrature squeezing, sub- 
and super-Poissonian photon statistics) of the constructed sp{N, R) algebra related CS 
and their superpositions are discussed and illustrated by several graphics. Our analysis 
shows that photon number oscillations are not necessary characteristics of nonclassicality 
of quantum states (neither are they sufficient P^). We note the main difference between 
squeezing in intelligent SS JI], H, ^ and in cat-type SS and construct a second kind mul- 
timode squeeze operator as a map from CS \d) to a set of cat-type multimode SS. In the 
appendix several statements of the main text are proved. 



2 The Barut-Girardello coherent states 



The property of canonical CS |a) [|I0[ to be eigenstates of photon number lowering operator 
a, a\a) = a\a) {a is complex number, [a, a^] = 1) was extended by Barut and Girardello 
to the case of Weyl lowering operator of su{l,l) algebra. Here we briefly review 
some of their properties. The defining equation is 

K_\z;k) = z\z;k), (1) 

where z is (complex) eigenvalue and k is Bargman index. Here, and in [II], we introduced 
k = — $ as a second label of the state and replaced the BG z with z/\/2. For discrete 
series D^^\k) the parameter k takes the values ±1/2, ±1, . . .. The Cartan-Weyl basis 
operators K± = Ki± iK2, K3 of sm(1, 1) obey the relations 

[K3,K±] = ±K±, [K^,K+] = 2Ks, (2) 

with the Casimir operator C2 = K^^ — {1/2)[K_K^ + K^K_] = k{k — 1). The expansion 
of these states over the orthonormal basis of eigenstates \k + n,k) of {K^ln + k,k) = 
{n + k)\n + k,k), n = 0,1,2,...) is 



^; k) = AfBG{\z\,k) ^ , \n + k,k)= ArBG{\zlk)\\z; k), (3) 

„=o ^n\r{2k + n) 



MBG{\z\,k) = [T{2k)/oF,{2k-\z\^)]'^ 



/2fe-i(2|z|) 
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where oFi{c; z) is the confluent hypergeometric function, Iu{z) is the modified Bessel 
function of the first kind, and T{z) is the gamma function [^. The above BG states 
\z; k) are normahzed to unity. Their scalar product is 



{k;z\z';k) = oFi{2k; z*z') 
and they resolve the unity (the identity operator), 

/ dn{z, k)\\z; k){k; z\\ = Ife, dfx{z, k) = -\z\^''''K2k-ii2\z\) d\ (5) 

where K^{x) is the modified Bessel function of the third kind. Note that \\z; k) = 
Mbg\z', k), while in [|| these nonnormalized CS were denoted as T{2k)~^^'^\z) (note also the 
misprint in [§: in formula for the measure function a{r) one should replace Ki^^i(2\/2r) 

by K2<i,+i{2\/2r) ^j). Owing to the above overcompleteness property any state |\E') can 
be correctly represented by the analytic function 



,F,{2k;\z\^)oFi{2k; 



(4) 



FBG{z,k;^) = {k,z*\^)/^^BG{\zlk) = {k,z*m, (6) 
which is of the growth (1,1). The orthonormalized states \k + n,k) are represented by 



monomials z"^ / n\T {2k + n) (we note a misprint in these monomials in 0] and 0: k 
should be replaced by 2k). The operators K± and act in the space Tik of analytic 
functions Fbg{z, k) as linear differential operators 

K+ = z, K^ = 2k^ + z^, Ks = k + z-^. (7) 
dz dz"^ dz 

This analytic representation has been used to explicitly construct eigenstates \z,u^v;k) 
of complex combinations uK_ + vKj^ in paper [|l|. 

Barut and Girardello have established their continuous representation for the discrete 
series D^{k), k = ±1/2, ±1,.... However, by inspection of their construction one can 
easily see that it also holds for reducible representations and forl/2>|/c|>0 - one only 
has to keep in mind that the quantity 1^ in the overcompleteness relation (^ is the identity 
operator in the subspace Hk, where su{l, 1) acts irreducihly. The proof consists of two 
observations (for concreteness we take D^{k)): (a) The expansions (|^) are convergent and 
represent normalized states for k > 0, provided \k,k + n) are orthonormalized; (b) The 
BG measure dfj,{z, k) resolves the unity operators by means of \z] k) for A; > provided 
the orthonormalized set of \k,k + n) is complete. 

It is well known that the s-u(l, 1) algebra has one- and two- mode quadratic boson 
representations, which are reducible in the spaces of states of one- and two-mode systems 
correspondingly. The one-mode realization of su(l, 1) is 

= \\ K2 = ^at2, Ks = ^{2a^a + 1). (8) 

Its quadratic Casimir operator C2 equals —3/16, C2 = — Kf — = k{k — 1), the 
Bargman index being k = 1/4, 3/4. The two-mode representation 

K_ = aia2, = alal, K3 = -{a\ai + 0202 + 1). (9) 
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is highly reducible (completely reducible), its irreducible components being just the rep- 
resentations from the discrete series D~^{k), k = 1/2, 1, . . .. The whole space Ti of the 
two-mode system states is a direct sum of the irreducible modules Tik- In these realiza- 
tions the operators uK_ + vK^, which were diagonalized in ref. |jl|, read ua"^ + f a^^ and 
uaia2 + va\a2- 

The Heisenberg-Weyl algebras hi and /12, spanned by 1, and 1, 

correspondingly, act irreducibly in the state spaces of one- and two-mode systems. The 
related families of CS \a) and |ai,a2) are overcomplete and realize the continuous repre- 
sentations, which proved to be very efficient |T^ . Therefore it is important to establish the 
relation between BG CS and the canonical CS representations. In the canonical CS rep- 
resentation every state l^') is represented by an entire analytic function -Fees (cn? 0^2; ^) 
of growth (1/2,2), 

Fbc5(ai,a2;^) = exp(i(|ai|2+ |a2r)) {a*i,a;\-^). (10) 

In the one-mode case Fccs{ci) = exp(i|Q;p)(a*|\E'). The eigenvalue properties of the 
BG CS and canonical CS and the realizations (||) and (|^) suggest that the canonical CS 
representation of a state |\E') G Tik should be obtained (up to a common factor) from 
its BG representation by means of substitution z = o? jl for the one- mode system and 
z = aia2 for the two-mode system, and this is the case. The corresponding relation 
between the two representations of the one- mode system states was written down in p|. 



Fccsioi; ^) = vr* 



^BGil(y^k = l) + ^aFBG{l(y^k = l] 



:iii 



If |\E') is even (odd) state, then the second (the first) term is vanishing. For the two- mode 
system states the relation between Fees and Fbg, defined above, is found in the form 
(proof in the appendix A.l) 



Fccsiai, «2; ^) = FBGiz,k = ^; ^) + ^ (af + af-') Fbg{z, k; ^), z = «i«2. (12) 

k>l 



Using these relations one can establish the coincidence between states, obtained in BG 
analytic representations and other familiar states. For example, the known one-mode 
even/odd CS \a)± coincide with the BG CS \z; 1/4) and \z; 3/4) while the gener- 

alized IS \z,u,v; k), constructed in using BG representation, for k = 1/4, 3/4 are the 
same as the eigenstates of ua^ + va^"^, constructed for real u, v in and for complex m, v 



in 0, |[ ^ 3 using the canonical CS representation. For A; > 1/2 the \z, u, v; k) with real 



u, V can be identified with two-mode SU{1, 1) states of |jT6|, All SU{1, 1) states 



of [32, 3^ can be found in the general family of su'"{l, 1) algebra related CS \z, u, v, w; k) 



constructed in 0, §|, 



In conclusion to this section it is worth noting that the SU(1, 1) group related CS |T0| 



provide another analytic (in the unit disk) |]2S[ representation of Hilbert space which has 
been shown p| to be related to the BG representation through a Laplace transform. It 
is also worth making a note about the notations: the BG CS are eigenstates of lowering 
operator = Ki — iK2, which belongs to the complexified algebra su'"{l, 1). Therefore 
we could denote such states as su'~"\l, 1) algebra related CS. However, usually when one 
deals with such simple complex combination as Weyl lowering/raising operators of an 
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algebra L {K± for sm(1, 1)) one writes L instead of (s-u(l, 1) instead of sm'"(1, 1)). For 
brevity we follow here this convention for BG CS for Lie algebras. Continuous families 
of eigenstates of general element of su^{l, 1) have been considered and called 1) 
algebraic CS or SU{1, 1) algebra eigenstates ||5|. An other motivation of the new term 
'algebra related CS' is the following property of the BG CS \z] k): unlike the h'^ algebra 
CS this family can not be represented in the form of group related CS either for the group 
SU{1, 1) or for the group of automorphysm Avit{su'^{l, I)) 3 SU{1, 1) [ pQ] . 



3 The BG CS for sp{N, R) 

The BG CS for semisimple Lie algebras can be naturally defined as eigenstates of mutually 
commuting Weyl lowering (or raising) operators Ea' (-EIO |]Ti1 ): 

Ea'\z) = Za'\z). (13) 

This definition can be extended to any algebra, where lowering/raising operators exist. 
We shall consider here the simple Lie algebra sp{N, R) (the symplectic algebra of rank 
and dimension N{2N + 1)). We redenote the Cartan-Weyl basis as Eij, Ejj, Hij = 
1,2, . . . , N, Eij = Eji, Ujj = Hji), and write the sp{N, R) commutation relations 

[E,„E,,] = [El,El]=0, 

[Eij, eIi] = 6jkHii + 6iiHjk + SikHji + 6jiHik, 

[Eij, Hki] = 6iiEjk + 6jiEik, (14) 

[El,Hki] = -6,kE]i-6,kEl 

[Hij, Hki] = 6iiHkj — 6jkHii. 



The BG CS for sp{N,R) are defined as eigenstates of E, 



Eij\{zki}) = Zi,\{zki}), t,j = 1,2,...,N. (15) 

Let us note that the Cartan subalgebra is spanned by Ha only and Hij^i are also Weyl 
lowering and raising operators as all Eij are: we have simply separated the mutually 
commuting lowering operators Eij. We shall construct explicitly the sp{N,R) BG CS for 
the quadratic boson rep, which is realized by means of the operators 

Eij = QiOj, Ejj = a|a], Hij = ^(a]ai + a^a]), (16) 

where N pairs of boson annihilation and creation operators. These operators act 

irreducibly in the subspaces H"^ spanned by the number states |ni, . . . , n^r) with even/odd 
ntot = ni + n2 + . . . + n^. The whole space H of the N mode system is a direct sum of 

The sp{N, C) is the complexification of sp{N, R) and therefore the Hermitian quadra- 
tures of the above operators span over C the sp{N, C) algebra. In the case of = 1 one 
obtains from ([T6|) the three operators K±^3 which close sp{l,R) ~ su{l, 1) (see equation 
(P)). We see that eigenstates of (the known even/odd states \a)± in quantum optics 



0) are sp{l, R) BG CS for k = 1/4, 3/4. 
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One general property of sp{N,R) CS for the representation (|T6D is that they 

depend effectively on complex parameters aj (not of A^^ + as one might expect). 
Indeed, using the boson commutation relations [oj, aj] = and the definition (|T5|) we 
easily derive 

^ij^kl Z'ik^jl ^il^jki iX^^ 

wherefrom we find the factorization of the eigenvalues 

Zij = aiOj, ai, aj G C. (18) 
Therefore in the above boson representation the definition (|15|) is rewritten as 

aiaj\{akai}) = aiaj\{akai}), i, j = 1,2, . . . , N. (19) 
The general solution to this system of equations is most easily obtained in the canonical 



CS representation |T0|. In Dirac notations the solution reads 



\{a,ai};C^,C^) = C4a)\a) + C4a)\ - a) = (20) 

where \d) are multimode canonical CS, a = (cti, a2, . . . , aAr) and C±{d) are arbitrary 
functions, subjected to the normalization condition (|ap = a- a= |aip + ... + |aArp) 



|C+(a)|2 + |C_(a)|2 + 2Re(C_Ci) iV(|a|) = 1, N{\d\) = {±d\ T a) = e-^l"!'. (21) 

Thus the families of states |a;C+,C_) represent the whole set of sp{N,R) BG CS for 
the representation (|l^). They have the form of macroscopic superpositions of multimode 
canonical CS. Macroscopic superposition of two canonical CS are also called Schrodinger 
cat states |ll4| , |15[, which we shall refer to as ordinary Schrodinger cat states. The set of 
( PPD is the most general family of superpositions of multimode CS \d) and | — c5). 



The large family of sp{N, R) CS ( pOD contains many known, in quantum optics, subsets 
of states 0, 15] and many other not yet studied. Let us point out some of the well 
known particular subsets of (pO]). The limiting cases of C_ = or C+ = recover 
the overcomplete family of multimode canonical CS, and and C_ = ±C+ produces the 
ordinary multimode even/odd CS |T5| . 

For the one-mode system (A^ = 1) several cases of the superpositions of two canonical 



CS ( PPD are thoroughly studied (for example, see for A^ = 1 and for any A^). 
Nevertheless, as far as we know, even in the one dimensional case no family of Schrodinger 
cat states was pointed out which is overcomplete in the strong sense in whole Ti. Here 
we provide such families for any A^. 
Consider in (|20|) the choice of 

C+ = cosip, C_ = isinyj, (22) 
which clearly satisfy the normcondition ( pT]) for any angle ip. 



\d-, ip) = cos ip\d) + i sin ip\ — a) . (23) 
In Fock basis (number states \ni, . . . ,12^)) we have the expansion 
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15; ^) = e-l"lV. ^ |n„ . . . , n.). (24) 

Using direct calculations we find that these states resolve the unity operator for any (p 
and thereby provide an analytic representation in the whole Ti, 

1 = —^ J d'^d\d;!f){(f;a\, d^a = dReaidlmai . . . dReaN dlmaN- (25) 
States |\E') are represented by functions 

/^(a,^) = el"^l'/2(^,a*|vl/), 
on which the operators aj and Oj act as 

d 

aj = P^aj, a] = -^v^^' (26) 

where acts as inversion operator with respect to ip: P^f{ip) = f\—ip). At <^ = 0, vr the 
multimode canonical CS representation ] = d/daj is recovered. 

The notations of ( pIS]) enable us to construct eigenstates of squared Weyl operators Ef^ 
(in any representation) as macroscopic superpositions of sp{N, R) BG CS in the form {zij 
are eigenvalues of Eij) 

\{z,i};D+,D^) = D4{z,,})\{zm}) + D4{z,,})\{-zm}), (27) 

where the functions D±{{zij}) have to be subjected to the normalization condition (sup- 
posing {{Zkl}\{Zkl}) = 1) 

+ + D^DI{{zm}\{-zm}) + D*iD+{{-zm}\{zm}) = 1. (28) 

In the quadratic boson representation ( [T6| ) these states take the form 



\{akai}; D+, D_) = D+\{a,a,}; C+, C_) + C+, C_) 

= \a-C+,C^,D+,D^) (29) 

and can be termed multimode squared amplitude Schrddinger cat states. They are expected 
to exhibit linear and quadratic squeezing and other nonclassical properties. In view of 
(pOD the states (^9]) are eventually expressed in terms of superpositions of four multimode 
canonical CS. 

In conclusion to this section we note that the overcomplete family of states \d]ip) 
admits n angles generalization: by means of n angles yj^, = 1, 2, . . . , n, n being positive 
integer, one can construct macroscopic superpositions of 2" CS \d) (or, equivalently, 
superpositions of 2"~^ sp{N, R CS of the type | a; </))), which are overcomplete and resolve 
the unity with respect to the same measure n~^d'^d, 

\d; (fii, . ..,(pn) = cosipn\d; v?!, • • • , V^n-i) +ismipn\ - d; ipi, . . . , V2„-l), (30) 
1 = ~?v / ^i,---, ^n){<^n, . . . , v^i; a|. (31) 
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In every component state in |c5; ipi, . . . , ipn) the parameters ctj are on a circle with radius 
\ai\. For n = we have CS |c5), for n = 1 the states (p3D are reproduced. |a; y^i, • • • , v^n) 
are easily seen to be eigenvectors of (ajaj)^" , and not of (ajOj)'", m < 2"~^, unless ipk 
are integer multiples of 7i/2. In the one mode case (A^ = 1) \a; ipi, . . . , ipn) are eigenstates 
of a^". Eigenstates of a^'' for A; = 1, 2, . . ., can be easily constructed as superpositions of 
sp{l,R) CS. Here we proved their overcompleteness for 2/c = 2" = 2,4,8,32.... Some 
eigenstates of powers of a'^, k > 2, have been considered in refs. [^]. Multicomponent 
macroscopic superpositions of canonical CS (one mode only so far) are intensively studied 
in quantum optics (with the final aim being the production of Fock states) ^ . 

The above result (31) is a particular case of a general theorem, proved in appendix A. 2, 
concerning the overcompleteness of common eigenstates of powers of non-Hermitian 
operators , j = I, . . . , N, n = 1,2, . . ., and valid for the case of A^ mode canonical CS 
and sp{N, R) BG CS as well. 



4 BG CS for the algebra u{p, q) 

The algebras u{p,q), p + q = N, are real forms of sl{N, C) and they are subalgebras of 
sp{N,R) [|ll|]. Therefore the BG CS for u{p,q) should be obtained from sp{N,R) CS by 
a suitable restrictions. In this section we consider these problems in greater detail in the 
boson representation ([T6|). 

The following subset of operators of (|l^) close the u{p, q) algebra (or m*" (p, q) if one 
consider non-Hermitian linear combinations of the operators below) [p!T[] , 

Eaf, = ao,a^, El^ = a^al, H^p = i(aj,a^ + a^al), H^^ = ^{a^a^ + a^a^), (32) 

where we adopted the notations a, P, j = 1, . . . ,p, fi, u = p + 1, . . . ,p + q, p + q = N 
(while k,l = 1,2, . . . , N). For p = 1 = q the three standard 5^(1, 1) operators K±, 
are K_ = E12 = aia2, = = ala^, = (a|ai + 0202 + l)/2. The subsets of 
Hermitian operators 

realize representations of compact subalgebras u{p) and u{q) correspondingly. The u{p, q) 



algebra (|32D acts irreducibly in the subspaces of eigenstates of the Hermitian operator L, 
L = T.Mi'2-i: = E - E H,, -{p- q)/2. (34) 

a ^ a ^ 

This is the linear-in-generators Casimir operator and the higher Casimirs here are 



expressed in terms of L |T3[ ) . Denoting the eigenvalue of L by / we have the expansion Ti = 
J2iZ-oo®'^i- The representations corresponding to ±/ are equivalent (but the subspaces 
T-L±i are orthogonal). We note that L = J2a ^a^ta ~ o)fj.o-tii and / = 0, ±1, . . .. 

The commuting Weyl lowering operators of u{p,q) are E^-^ = a^a^, 7 = 1,2,. . . ,p, 
fi = p + l,p + 2, . . . ,p + q = N. We have proved in the above that eigenvalues of the 
product of two boson destruction operators are factorized. Therefore the u{p, q) BG CS 
in the above boson representation can be defined as |{a/3ai/}; l,p, q). 
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a^a^\{af3au};l,p,q) = a^a^\{ai3au};l,p,q), 

(35) 

7 = 1,. iJ = p+l,...,p + q, 

where and are arbitrary complex numbers. We put \\d;l,p,q) = g), 
denoting by a nonnormalized (but normalizable) state, while is normalized to 
unity. Solutions to the above equations can be written in the form 

ni n]v_i "p-nq-i 

\\a;l,p,q)= 2^ I =\ni,...,nN-i;np-ng-l}, (36) 

hj,-n,=i ^nil... nN-il{np -n'^-l)\ 

where aj, i = 1,...,N, are arbitrary complex parameters, hp = Y.a'^ai = I]^n^ 



h' = hq — rij^ and I = hp — hg. In (|36|) summation is over all rii = 0, 1, 2, . . . provided 



Up — Uq = I = const. 



If we multiply q) by exp(— |ap/2) and sum over / we evidently get the normal- 

ized multimode CS \a) (for any pair p, g), 

oo 

\S,) = e-h\^? Y: \\a;l,p,q). (37) 

/=— oo 

The last equality suggests that the states ||c5; l,p, q) form overcomplete families in Tii for 
every p, q. This is the case: using the overcompleteness of formula (|37| ) and the 
orthogonality relations 

(p, g, g) = for I'^l, (38) 
one obtains the resolution of unity in Tii in terms of the m(p, g) CS g), 

/" (i/i(a)||a; g)(p, g, /; a|| = 1;, dfi{d) = — rTe~'"'^(i^a. (39) 



Now we note that in u{p, q) CS ( PB| ) one complex parameter, say oat, can be absorbed 
into the normalization factor by redefining the rest as 

zi = aiaN, . . .,Zp = apUN, ^p+i = ftp+i/ "tv, • • • , zn-i = aN-i/aN- (40) 
Then we can write ||c5; l,p, q) = a~M\z] l,p, q) and 

Zi^ z^~^ 

\%l,P,q)= I , ^ - =\ni,...,nN-uhp-h'q-l), (41) 

where z = {zi, . . . , zn-i)- The states \\z; l,p, q) are normalizable in view of 

oo 



/=— oo 



which stems from ( |37D and ( |38D . The normalized states | /, p, g) are | Z, p, g) = A/" || /, p, g) , 
A/" being the normalization constant. 
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The family {\\z; l,p, q)} is overcomplete in Tii and the resolution of unity reads {d^z = 
nf dReZi dlmzi = {onI^'^'^-^'P^ Uf'^ dReai dlmai) 

h = J dn{z;l,p,q)\\z;l,p,q){q,p,l;z\\, 

^ ^ , (42) 

dfi{z,l,p,q) = F{\zp\, \zg\;l,p,q)d^z, 

where |5pp = {zi]"^ + . . . + {zpl"^, = l^^p+ip + . . . + I^at-iP, the measure weight 

function being 



F{\^pl\~Zg\;l,p,q) = — I d^aNlaNl'^'^''-"''^ 



exp 



I pi I l£^|2u |2 I U |2 



+ \Z„\ Oat + Oat 



2 ' 



(43) 

One can prove that the above measure is unique in the class of smooth functions of 
|zi|, . . . , |-2Ar_i| (see appendix A. 3). Thus the explicit form of u{p, q) BG CS is 

\z]l,p,q) = A^d^;!!,. . . \zN-i\;l,P,q)\\z;l,p,q), (44) 



where \\z; l,p, q) take the form of superposition (|41|) of multimode Fock states with fixed 
value / of the difference number operator L, equation (p^). 

Let us note some known particular cases of the u{p, q) BG CS (^TJ). Recently the case 
of g = 1 and negative Z, — / > (then p = — 1, = 0, = and Zp = z) has been 
considered by Fujii and Funahashi [0]. Their resolution unity measure (in Tii) reads 

d^^\z) = F'{\zll,p^)d''z, F' = ^^4— ^"'-P+i(2|^), (45) 

where Ky{z) is the modified Bessel function of the third kind |^. F'{\z\,l,p,l) and 
F(|z|,/,p, 1) do not depend on phases of Zi and are smooth functions of . . . , 
i.e. all order derivatives are finite. In appendix A. 3 we prove that the resolution unity 
measures for u{p,q) CS are unique within such a class of functions, i.e. F'{\z\,l,p,l) 
and -F(|i1, l,p, 1) should coincide. Then using the analyticity property of Bessel functions 
K^(z) ||2^ we establish (proof in appendix A. 4) the following integral representation for 
K„{z) with 1^ = 0, ±1, . . . and Kez > 0, 

K^{2z) = -z-" r dxx''-^e-^''+''^''\ (46) 
2 Jo 

For p = 1, q = 1 our states \z] l,p, q) recover (as the states of 0] do) the BG CS \z] k) for 
the series D~^{k) of su{l, 1) the Bargman index k being expressed in terms oi I as k = 
(1 + |/|)/2. The irreps with ±1 are equivalent, however the states \z; /, 1, 1) and \z; —1, 1, 1) 
are different as one can see from their definition (l4l|) (moreover, they are orthogonal). 
Thus our states \z; ±1, 1, 1) represent two equivalent but different realizations of BG CS 
\z; k) for = (1 + |/|)/2 = 1/2, 1, . . .. The exact identification is \z] / < 0, 1, 1) = \z] k), 
||2;/>0,1,1) = z'^''-^\\z;k). 

The pair CS in quantum optics \C,q) [0] (defined as eigenstates of aia2 with a\ai — 
(^2(^2 = Q = const.) appear as u(l, 1) BG CS \z] k) in the two-mode representation (A^ = 2 
in (|32|) ). The identifications is |C, q) = \z] /, 1, 1), i.e. the Agarwal C and q are equal to our 
z and / correspondingly. In view of Eqs. (^) - (|i3| ) the pair CS are overcomplete in the 
subspaces Tii. Our \z; l,p, q) can be regarded as a generalization of |C, q) to the A^-mode 
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boson system: |z; /, p, g) are invariant under the annihilation of pairs of two different mode 
bosons, one from the first p modes, and the other from the last q modes. Note that in the 
sp{N,R) CS |a, C_,C+) there is no such restriction - these are the most general states, 
which are invariant under the annihilation of any pair of bosons. 

The sp{N, R) CS |a, C_, C+) can be decomposed in terms of u{p, q) CS \z] l,p, q) with 
different /. For N = 2 this decomposition reads 

oo oo 

\a; C_, C+) = a[Q\z; 1, 1) + ^ /, 1, 1), (47) 

1=0 1=1 

Ci = C+ + z = aiaa- 

In analogy to the case of sp{N,R), considered in section 3 we introduce the u(p,q) 
multimode squared amplitude cat states \z] l,p, g; D^, DJ) as macroscopic superpositions 
of g) BG-type CS \z;l,p,q), 

\z;l,P,q;D+,D^) = D+\z;l,p,q) + D^ \ - z;l,p,q), (48) 

which are expected to exhibit squared amplitude squeezing and other nonclassical prop- 
erties. In the particular cases of = 2, D_ = exp{i(f)) the states (^Sf ) recover the 
two-mode Schrodinger cat states, considered recently in pO[ . 



5 Statistical properties of the mode sp{N, R) BG 
CS and their macroscopic superpositions 

In the present section we consider some general statistical properties of the constructed 
sp{N, R) algebra related CS and their superpositions and discuss in greater detail some 
new subsets of this large family. 

All sp{N, R) BG type CS minimize the Robertson multidimensional uncertainty re- 
lation ||38| for the Hermitian quadratures Xij, Yij of mutually commuting Weyl lowering 
operators E^, since they are eigenstates of all Eij (Proposition 3 of ref. 0), 

det a({X,„ Y,,}; a, C_, C+) = det C({X,„ a, C_, C+), (49) 

where a is the matrix of second moments of all observables Xij,Yij (the uncertainty 
matrix) and C is the antisymmetric matrix of all mean commutators of Xij,Yij times 
{—i/2). The number of commuting Eij is equal to (A^^ + N)/2. Robertson inequality for 
n observables Xj, j = 1,2, . . . ,n, reads 

det ct{{Xj}; ^) > det C{{Xj}; ^), (50) 

and for a pair of two operators X, Y it reduces to the Schrodinger one, A^XA'^Y — aj^y — 
|([X,F])|V4, where axx = {XY + YX)/2 - {X){Y) (for greater detail see for example, 
m, 0). In all sp{N, R) BG CS the covariances of Xij and Yij are vanishing, but those of 
Xij and Xki are not, i.e. the matrix cr({Xjj, Yij}; a, C_, C+) is not diagonal. 

The subset of Schrodinger cats \(y,ip), equation (p3D, possess several remarkable prop- 
erties: 
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(a) They are overcomplete in the whole Hilbert space (see equation (p5D); 

(b) The photon statistics in every mode is Poissonian for any a and if. This fol- 
lows immediately from the expansion (^) in terms of multimode number states \n) = 
\ni, . . . ,nAr); 

(c) the states \d;ip) can exhibit squeezing in the quadratures pj, qj (for example, for 
\d\ close/equal to \ai\ = 0.5, = 7r/4 and arga, around n7r/2, n = 0, 1, . . ., the minimal 
value of Api and Aqi being equal to 0.316 - see the graphics /i on Fig. 1 ); 

(d) these states are physically coherent ('true coherent') since they satisfy the condition 
of full second-order coherence of the field The latter property again follows from 



equation (p^, which is of the form of generalized CS of Glauber and Titulaer [p9[] . 

This interesting subfamily {la;^?)} of sp{N,R) BG CS can be generated from the 
familiar multimode canonical CS by means of the following operator 

S{^) = exp (2(-l)>) : = S{v)\a). (51) 

where h = a\ai + . . . a}j^aN is the total number operator. As strange as it may seem S{(p) 
is well defined for any angle (p and is unitary. On any state |\E') its action is 

5(<^)|vI/) = e*'^||vI/), + e-^^||vl/)„, (52) 

where ||\l')e,o are the projections of |\E') on even/odd subspaces The operator (— 1)" 
is Hermitian and {—l)^ip may be regarded as a sort of nonlinear multimode interaction. 

In the classification scheme of [|TB] the (one-mode) states which possess the above 
properties (b) and (c) fall into the subclass of the weakly nonclassical states. In this scheme 
the nonclassical states are subdivided into weakly nonclassical and strongly nonclassical 
depending on the pointwise nonnegativity or nonpositivity of the phase averaged 'P{I) 
Glauber-Sudarshan diagonal representation P{/3), I = /3 = y/l exp{ii!}), 

V{I) = P(re*'^)d¥, r = Vl = \p\. (53) 



2n Jo 

If V{I) < for some values of / the state is strongly nonclassical (then also -P(/3) < 
for some values of (3) and if V{I) > but -P(/5) ^ the state is said to be weakly 
nonclassical |T3]. The set of classical states (i.e. P{l3) > 0) is not subdivided. Criteria 



for phase-insensitive nonclassicality of single-mode states were also studied in |n 

The family of sp{N,R) BG CS equation (p3D, consists of classical (at ip = 

0, ±7r/2, tt) and weakly nonclassical states (for ip ^ 0, ivr, vr) for every mode since the mul- 
timode photon distribution in these states is a product of one-mode Poisson distributions. 
There are not strongly nonclassical states in this family. Note that at ip = 0, ±vr/2, vr the 
states |a; (p) are the CS \a) or | — «), and at = 7r/4 they coincide with the Yurke-Stoler 
states [Q. The states with Gaussian Wigner function are either classical or strongly 
nonclassical [T^ and strongly nonclassical states from the latter family all have posi- 



tive Mandel Q factor (super-Poissonian statistics) ^ [Q = (A^ra — (n))/(n), where 
n = a^a]. 

Along these lines we note that in the family of weakly nonclassical states \a; p>) there 
are states which exhibit quadrature squeezing (graphics /i on figure 1). Conversely, there 
exist strongly nonclassical states (for example, in the family ■?/'), defined below) 
which do not exhibit squeezing of the quadratures of either a or (nor the Q factor is 
negative). Moreover, among the one- mode \a, 0, ip) there are states with Q = which are 
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squeezed or not squeezed, but their photon statistics is not Poissonian (see graphics on 
figures 4 and 5). These examples show that Q = is not a sufficient condition either for 
a statistics to be Poissonian or for a state to be classicaL 

The quadrature squeezing and/or Q < are sufficient conditions for the nonclassicahty 
of the corresponding states [41]. However, they are neither sufficient nor necessary for the 
strong nonclassicahty as demonstrated below. 

In |]18[ a simple sufficient condition for strong nonclassicahty of the states (i.e. for 
non positivity of the phase smeared diagonal P representation V{I)) is given in terms of 
photon number distributions Pn, 

In '■= {n + l)pn~iPn+i — < fo^ somc w > 0. (54) 



The distribution p„ is expressed in terms of V{I) as JTS 



Pn 



POO 1 

/ d/P(/K(^"^)(J), pJ^-)(/) = re-^. (55) 
Jo nl 

Distributions pn which can be represented in the above form with V > {V ^ 0) were 
recently defined as classical (nonclassical) . Nonclassicahty of p„ means strong non- 



classicahty of the corresponding states. 

Among sp{N, R) BG CS there are strongly nonclassical states as well (the definition 
of strong nonclassicahty for multimode states is discussed below), such as, for example, 
the cat states 



\a,(j)) =Ar{\a) + e'^\-a)), (56) 
the normalization constant being 

Af= (2(l + cos0e-2*='))"^ =Ar(f,0). 

In the case of iV = 1 the states ( |5^ ) have been discussed, for example, in [^, ^ and in 
the fifth paper of [|1^. The probability for totally n photons in \d, (p) (irrespective which 
mode they belong to) , n = rii + ... + tin , is found as 



p„(f,0)=iV(f,</))V^'— s„(0), s„(0) = 2(l + (-l)"cos0), f=\al (57) 

nl 

and the function Z„(f, 0) takes the form 

L{r, 0) = ln{r, (p) (s„_i(0)s„+i(0) - s^(0)) , (58) 

where the function 

/„ = iV(f,0)V^-— -, 

nl[n — 1)1 

is nonnegative. The nonnegative factor s„(0) is seen to be a bounded and oscillating 
function of both and n, Sn{4>) = Sn+2{4>)- Then for every ^ ±7r/2 the combination 
s„_i(0)s„+i(0) — s^(0) is negative for all n for which (— l)"cos0 > 0. Noting that the 
total photon number distribution p„(f, 0), equation (|57D, coincides with that for the one- 
mode states |a,0), \a\ = r, we conclude that all one-mode states |a,0 7^ ±7r/2) are 
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strongly nonclassical (the strong nonclassicality of the one-mode states |a, 0) was proved 
also in the very recent E-print |2^). 



One way of generalizing the notion of strong nonclassicality to multimode states is to 
apply the above definition to the total photon number distribution, the other one is to 
require this for every mode. One can easily verify, that in the N mode states |a, 0) the 
conditional photon distributions Pni,...,ni,...njv('5) 0) for individual mode i {uk^i being fixed) 
also obey the inequality (|5^). Thus |a, 7^ ='=7r/2) are strongly nonclassical according to 
both criteria. 

Consider now the squared amplitude quadrature squeezing [O in the multimode states. 



We first note that the BG-type CS for any Lie algebra cannot exhibit squeezing of the 
quadratures and Yij of Weyl operators Eij since here the variances of Xij and Yij 
are equal which stems from their eigenvalue property ([T5|) [|I]. In the quadratic boson 
representation Eij = a^aj and X^j (or Yij) squeezing is multimode squared amplitude 
squeezing. The quadrature X^j {Yij) of ajOj is said squeezed in a state if the variance 
AXjj (AFjj) is less than its value in the ground state |0). Thus quadratic field squeezing 
does not occur in |a;C_,C+). We shall see that macroscopic superpositions of two such 
states do exhibit quadratic squeezing. But let us first make some general remarks about 
SS of two and several observables. 

Squeezing of the two quadratures X and F of a non-Hermitian operator A (for defi- 
niteness we write A = X + lY) can be achieved in two ways: 

(a) in the eigenstates |z, u, v) of complex combination uA+vA^ ( generalized intelligent 
states) 0; 

(b) in the eigenstates \zY'^^ of A^ (generalized cat states). 

The first possibility was proved and demonstrated (on the examples of SU{1, 1) and 
SU{2) generators in series D~^{k) and D{j)) in [0, ^. These SS minimize the Schrodinger 
inequality and therefore were called Schrodinger (or generalized) IS. A particular case 
of SS of type (a) are the SS for general systems introduced as states minimizing 
the Heisenberg inequality, which is a particular case of that of Schrodinger. The second 
possibility (b) can be easily proved by calculations using the eigenvalue condition of A"^ 
and taking into account the Schrodinger relation. This can also be checked directly on 
the example of the following two type of superposition states 

\z; (f) = cos (p\z) + i sin (p\ — z), \z, (f)) = Mdz) + e'^'^\ — z)), (59) 

where \z) are eigenstates of A, A\ ± z) = ±z\ ± z). These states can exhibit squeezing 
according to the stronger criterion, given in [|I| (see also below). \z:,(p) and \z,(j)) are 
eigenstates of (and not of A, unless = rni /2, = 0, 1, . . .). Eigenstates l-z)*-^-* of A 
which are not eigenstates of A are superpositions of | ± z). Therefore SS of type (b) are 
cat states. \z\ (p) and \z, 0) in (|59|) are examples of such SS for any A for which eigenstates 
I ± do exist. In fact first- and higher order squeezing of X and Y can occur in states 
which are eigenvectors of A^ for any n>l and such eigenvectors can be easily expressed 
as discrete superpositions of several \z). 

The operator S{u,v) which transform the nonsqueezed \z) to the SS of type (a), 
\z, u, v), was defined in |^ as generalized squeeze operator {if A = a then S{u, v) is the 
known canonical squeeze operator |^). Having established that eigenstates l^;)*-^-* of 



A"^ can universally exhibit squeezing of the quadrature of A we can define, in analogy to 
the previous case, a second kind of squeeze operator Sjj by means of the relation 
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\z)^^^ =Sn\z). (60) 

We can point out an example of such a squeeze operator - that is the operator S{(p) of 
equation (^TJ). It maps the multimode CS \a) to the weakly nonclassical states {sp{N, R) 
BG CS) 1 5; Lp) which are eigenstates of a^aj (the sp{N, R) BG CS) and do exhibit quadra- 
ture squeezing (see the graphics /i on figure 1). 

The main difference between the above two types of SS is the following: SS of type (a) 
can exhibit arbitrary strong squeezing of X or Y, while the squeezing in SS of type (b) 
is always bounded, since eigenstates of ~ (X + iY)'^ can never tend to an eigenstate 
of X or Y. Family of states in which arbitrary strong squeezing ('ideal squeezing') of 
X or F is possible could be called ideal X-Y SS. So the Schrodinger IS, in particular 
the canonical SS are ideal p-q SS. We follow the definition of X-Y SS according 
to ref. a state |^) is X-Y SS if AX < Aq or Ay < Aq, where Aq is the lowest 
level at which the equality AX = AY can be maintained. The lowest level is reached on 
some eigenstate \zq) of A. For the quadratures of a^, k=l,2,. . . , and {aiajY the lowest 
level is reached in the ground state |0). Linear and/or quadratic squeezing in ideal one- 
mode squared amplitude SS (eigenstates \z;u,v) of ua? + va?) is considered in papers 
(for different ranges of parameters u, v) |T^, |3^, ||, The diagonalization of 

ua'^ + va^'^ for k > 2 is discussed in the very recent E-print |^ . 

A family of states in which squeezing of quadratures of any product ajaj, i,j = 
1, 2, . . . , X, can occur should be called a family of multimode squared amplitude SS. Exam- 
ple of such multimode SS is given by the Robertson IS 0, which should be eigenstates of 
complex combinations Uki-ijaiaj + Vki^ijala^j (summation over repeated indices). These are 
ideal multimode squared amplitude SS. Multimode quadratic SS of type (b) are defined 
as eigenstates of all squared products {aiaj)^. They take the form (^9]) . 

Next we consider cat type squared amplitude SS. One example of two-mode cat-type 
second-order SS is considered in paper |^|, which, however, was examined for ordinary 



squeezing only. Here we provide examples of multimode cat-type SS which can exhibit 
both quadratic and linear squeezing and other interesting statistical properties. Such SS 
are the following macroscopic superpositions |a, 0, ip) of the sp{N, R) algebraic CS \a, 0) 
(|a, 0) are defined in equation (^)): 

|a,0,V^) = Ar(|a,0) + e'^| -a,0)), (61) 
where J\f is the normalization constant, which obey (|21| ) and has the form 

Af{r,^,ij) = (l + 2A/"2e-^'' (cos0 cos(f2 - + V^) + cos(f2 + - ^))) , (62) 



Af being given in equation ( p^ ) and r = \d\ = yjaiP -|- |a;2p + . . . -|- |aArP. 

We demonstrate the quadrature squeezing on the example of individual mode op- 
erators tti (linear squeezing) and aj (quadratic squeezing). Note that 1 5, 0,-0) are not 
factorized over the different modes. The variances Api and Aqi of the quadratures of i 
mode annihilation operator a^, = (gj + ipi)/\^ are 

A^Pi{f, n, 9i, 0, 7/-) = ^ + {aUi) - Re(a^) - 2{lm{a,)f, 
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A^iir, n, 9i, cj>,^) = - + (aUi) + Re{a^,) - 2{Re{ai)y, (63) 
where Vi = \ai\, 9i = argaj and 

(fli) = -2Q^iM'^Af^e-^' sincf) {1 + i) (e"^" + cos(f2_0+^) + sin(f2_0+^)) , 
(ajai) = 4r^N"^J\f^ (l - cos0e"^^' - e"^' (cos0 sin(f^-0+'0) + sin(f^+0-'0))) . (64) 

As functions of 9i the variances of pi and qi oscillate with period tt and Api{9i + 7r/2) = 
Aqi{9i). Linear squeezing is exhibited in states, for example, |a, O,^') with = 0.05, f 
close to Tj, 9i = 7r/4, (f) — and -0 around 3.152 (see figure 1). Maximal Pi {qi) squeezing 
is obtained when r = Vi (i.e. when only one mode is excited). Here A^pj > 0.275 = 
A2pi(0.05, 0.05, 7r/4, 0,3.131) = A2gi(0.05, 0.05, 7r/4, 0, 3.153). When f is increasing the 
graphics of Api and Ag^ (as functions of the angles (p and ip) become smoother and tend 
to a constant value, independent of the superposition parameters (p and ip. In the above 
states \a,0,ilj) the Mandel factor Qi for the mode i is negative in the vicinity oi ip = n 
only. By its definition the quantity coincides with the intensity of the field (the total 
mean number of photons (a^a) = I]f^(ajaj)) in the multimode CS \d). The field intensity 
in the multimode superposition states \d,(j),ijj) reads 

(a^a) = Af'^Af'^N"^ (l-cos0e"^^'-e"^' (cos0 sin(f^-0+V') + sin(f^ + 0-V'))) . (65) 

We see that (a^a) is an increasing function of f, f = 

The variances of the quadratures Xi, 1^ of the squared individual mode amplitude af , 
a? = {Xi + iYi)/y/2, in |a, 0, -0) is easily obtained in the form 

A'Xi{r, n, 9i, 0, V) = 1 + 2(aja,) + (af a,^) + rf cos(4^,) - 2(Re(a,^))^ 
AV,(f, n, 9i, 0, V') = 1 + 2(aja,) + (af a^) - rf sin(4^,) - 2(Im(a,2))2. (66) 

where 

(a^) = -Ua\N'^M'^e~''^ (cos0sin(f^ - + -0) - sin(f^ + - -0)) , 
(of ) = 2r\}P (l - 2A'^e-*^' (cos cos(f2 - + V) + cos(f2 + - i/-))) . (67) 

As functions of the angle 9i the variances of Xi and Yi oscillate with period 7r/2 and 
AXi{9i^T^I^) = AYi{9i). 
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Fig. 1. Amplitude quadrature squeezing in 
the sp{N, R) BG CS |a; (p) and the superpo- 
sitions \d,(f),ilj), equation (61), for = |ai| = 0.05. 

Aqi = Aqi{f,ri,ei,(p,tjj), f=\a\, qi = qi or pi. 
fi = A'^qi{ri,ri,-^,ip), f2= A'^pi{ri,ri, 1,0,^;), 
f3 = A\iri, n, f , 0, V), h = A^pi{iri, n, f , 0, V), 
I a, (f) are weakly nonclassical states for eve- 
ry mode, \d, (f), ip) are strongly nonclassical. 



Fig. 2. Squared amplitude quadrature squeezing 
in superpositions states |a, ^, ■0) for = 0.8. 

AXi = AXi(f,r„ei,0,V), Xf = Xi or Yi, r = \d\. 
gi = A^X,{ri,ri,l,0,i;) = A^Yi{ri,ri,-^,0,ij), 
52 = A2x,(1, r„ f , 0, V), 53 = A2X,(1.2, r„ f , 0, V), 
g^ = 2A^Pi{ri, n, f , 0, i;) = 2A^qi{ri, n, -f , 0, V). 
Joint X and p (or Y and qf) squeezing occurs in 
the interval 6.4 < ^ < 7.4. 
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Fig. 3. Probabilities pn{f,(f),^j) to 
find n photons in the multimode super- 
position states \d, (p, Tp) as functions of 
f = \d\ for different values of and 
V'- Po = Po(r, f,7r), pi = pi(f,7r,-f), 
P2 = P2(r,0,7r), P3 = P3('',7r, f), P4 = 

P4(r, f , f ), P5 = Pb{f, TT, -f ). 



Fig. 4. Oscillating photon number distri- 
butions Pnir-, <pj "P) ill strongly nonclassical 
states \d, (p, ip) for different values of f, (f) 
and ip. pnl = Pn(0.8,0, 7.3) (Q > 0, Ap > 
0.38, AX > 0.73), p„2 = p„(2.2, vr, -f ) 
(Q < 0), p„3 = pn(0.55, 5.0914,0) (Q = 
-0, Ap > 0.38, AX > 1). 



The variances AXi and AYi are squeezed if they are less than their value of 1 in 
the ground state |0). This holds, for example, in states \d,0,ilj) with f close/equal to 
^ 1; 01 = mr/A, = and ijj around zero, the minimal value of A^Xj and A'^Yi (at 
f — Ti — 0.88) being equal to 0.69. In the states |0.8e''=''^/^, 0, ■V^) linear and quadratic 
squeezing can occur simultaneously (see graphics gi and 514 on figure 2: joint X^- and pi- 
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{Yi- and qi-) squeezing occurs in the interval 6.4 <■?/'< 7.4). In the above interval Qi > 0, 
where Qi is the Mandel factor for the individual mode i. 

On figure 2 graphics are shown of A^Xj(f, rj, ir/A, 0, ip) as a function of the angle ip for 
fixed Tj = 0.8, and three different values of the total excitation parameter r, r = 0.8 = 
(i.e. only mode i excited, graphics gi), r = 1 (graphics (72) and r = 1.2 (graphics g^). 
One sees, that graphics of AXj(-?/') become rapidly smoother and tend to a constant value 
when f is increasing. 

An important statistical property of all states \d,(j),tp) is that they are strongly non- 
classical in the sense of definition of |[T^ (discussed above), which we apply here to the 
total photon number (and to the conditional individual mode number) distribution in 
the multimode states. The total photon number distribution p„(r, 0, ip) takes the form 
similar to that of equation (|57|) , 



1 + (-l)"e^* + i"'e'^ + (-ire'^^-'^'> " 



(6^ 



where 



■2/';; 



r, 0, ip)Af' 



2^, 



The factor Sn{(f),'ip) is bounded from above and as a function on n oscillates with period 
4. Therefore the inequality (Q) is satisfied in all states |a, 0, ■0) for those n for which Sn 
reaches its local maximum. This proves that all |a, 0, -0) are strongly nonclassical. Note 
that the factors for conditional distribution of Ui {uk^i fixed) also satisfy the inequality 
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0.3 
0.2 
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p„4 =Poisson ■)<■■■ 
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Fig. 5. Nonoscillating photon number 
distributions in strongly nonclassical states 
I a, 0,0) for different values of f = \a\, 
and V- Pnl = Pn(0.55, 2.246,0) 
(Q < 0, Aq > 0.5, AX > 1), p„2 = 
p„(0.55, 2.33,0) (Q < 0, Ag > 0.43, AX > 
1), pn3 = p„(0.55, 2.234384,0) (Q = +0, 
Aq > 0.5, AX > 1), p„4 = Poisson distri- 
bution with (a^a) = 0.685 as in p„3. 



n 



As in the case of |a, 0) here again p„(r, 0, tp) coincides with the probability to find n 
photons in the one-mode states |a,0, -0), |a| = r, a = fe*^. The distributions pn{r,4>,4') 
does not depend on 6. It can be oscillating or nonoscillating and with positive, negative 
or vanishing individual mode Q factor. No definite relations exist between the sign of 
Q, the photon number oscillations and the amplitude quadrature squeezing: all possi- 
ble combinations of these three properties can be found in strongly nonclassical states 
|a, 0, -0). In figures 4 and 5 representative graphics of oscillating (figure 4) and nonoscil- 
lating (figure 5) photon distributions are shown. The sign of the corresponding Q and the 
inequalities for Aq{6) and AX (6) for each graphics are also given. In the recent E-print 
p4| examples of classical states with oscillating photon distributions were pointed out. 
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Thus photon number oscillations are neither necessary nor sufficient for nonclassicality of 
quantum states. 

The Q factor is bounded, Q > — 1, and when Q = —1 then the variance An of h is 
vanishing. This means that the corresponding state is an eigenstate \n) of n (a Fock 
state) and Pn(|^)) = 1- In figure 3 photon probabilities Pn{f,(j),ip), n = 0,1,2,3, are 
shown as functions of f for several values of and ip. At f — one obtains Pn = ^- For 
> 1 this yields the finite superpositions of multimode Fock states |n) , no + . . . + = n, 
and for the one mode case, = 1, - the number state |n) with n = 0,l,2orn = 3. We 
see from figure 3 that practically the states |a, 0, ip) with the corresponding 0, ip coincide 
with Fock states |n), n = 1, 2, 3, for |a| < 0.5 (then p„ > 0.99995). In multimode case for 
|a| < 0.5 the specific form of superposition of several \n) depends on the specific values of 
+ . . . + \aM\ < 0.5. If ak^i = then all n photons/bosons are in the mode i, i.e. 
the Fock state is |n) = |0, . . . , n^ = n, 0, . . . , 0). There is a growing interest in obtaining 
Fock states from macroscopic superpositions of (so far mainly one mode) CS \a) (see |^ 



and references therein). Here we provided example, probably the first one, of obtaining 
Fock states of multimode systems. 



6 Concluding remarks 

We have constructed and discussed some properties of sp{N, R) and u{p, q) algebraic (al- 
gebra related) CS in the quadratic boson representation. These states are a generalization 
of the sm(1, 1) CS of Barut and Girardello and are constructed as eigenstates of all mu- 
tually commuting Weyl lowering operators. The quadratic boson realizations of sp{N, R) 
and u{p, q) are reducible. Therefore the corresponding group related CS are not over- 
complete in the whole Hilbert space of states H. The BG-type CS are very large sets 
and afford the possibility to resolve the unity operator in Ti. by means of some subsets. 
We pointed out such subsets of the sp{N, R) algebra related CS (and their superpositions 
as well) and wrote down the relations between the established u{p, q) CS representations 
and the familiar A^ mode canonical CS representation, in particular between the su(l, 1) 
BG CS and the two-mode canonical CS representations. 

The new states can exhibit interesting statistical properties, such as amplitude quadra- 
ture squeezing, sub- and super-Poissonian photon statistics and oscillations in photon 
number distributions. All states from the overcomplete subfamily |a;v?) of the sp{N,R) 
BG-type CS are weakly nonclassical and (some of them) can exhibit amplitude 



quadrature squeezing as well. Strongly nonclassical [0 sp{N, R) algebra related CS 
were also pointed out. 

Noting that the BG type CS \z) can not exhibit squeezing of the quadratures of Weyl 
generators Eij we anticipated that such squeezing should occur in eigenstates of E^, 
m > 2, which for Efj = {aiajY are called multimode squared amplitude Schrodinger 
cat states. Squared amplitude squeezing in the individual modes is demonstrated in the 
superpositions |a, 0, ip) of two sp{N, R) CS. These are strongly nonclassical states and at 
small \d\ {\d\ < 0.5) and for specific values of the angles 0, ip practically coincide with 
superpositions of several multimode Fock states |n) with total number of photons/bosons 
n = 1, 2 or n = 3. If a^^j = then all n photons are of the mode i, i.e. we have a 
single multimode Fock state. The Fock state engineering via discrete superpositions of 
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canonical CS \a) is of current interest in the literature (one mode mainly). We have 
shown that discrete superpositions of multimode canonical CS are naturally encompassed 
in the framework of sp{n, R) BG-type CS and their linear combinations. The weakly 
nonclassical sp{N, R) BG CS \d; ip) can be generated from CS |a) by means of the second 
kind (unitary) squeeze operator. In greater detail this should be considered elsewhere. 

Acknowledgment. The work was partially supported by Bulgarian Science Founda- 
tion under contract no F-559. 



7 Appendix 

7.1 A.l. The correspondence rule between the mode canon- 
ical CS and u{p,q) BG-type CS representations 

The multimode CS \d) are overcomplete in the mode boson system Hilbert space 7i, 
spanned by the number states \n) = \ni, . . . ,ni\i). In canonical CS representation a state 
|\&) is represented by analytic function Fccs{c(] ^) of ^ variables ai, i = 1, . . . , N, 

Fccs{a,^) = {a'^\m, ||«) = E " V , 1^) (69) 

The u{p,q) BG-type CS \\z]l,p,q), equation (|4l|) , are overcomplete in subspaces Hi, 
satisfying the resolution unity equation (^2]). Hereafter if a state l"^) E Hi C H, then in 
the u{p, q) CS representation this state is represented by the analytic functions of — 1 
variables Zk, k = 1, . . . , N — 1, 

Fi{p,q,z--^) = {q,p,l-z*\\'^). (70) 

The relation between the two representatives of \^) immediately follows from the expan- 
sion (|37|) and equations (|40|) and (0): 

oo 

Fccs(«, ^)= E Oi~NPi{P^1^^'^'^)^ Zfc given by equation (40). (71) 

l=—oo 

This formula is efficient for the transition from {Fi} to Fqcs if one knows the represen- 
tatives Fi{p,q,z; \E'). In the opposite direction the transition formula is easily obtained 
from ([TOD, { ^7\j and the orthogonality between Hi and Hi'^i, 

1 r 

Fi{p,q,z';^) = ^ j cfaa'^{q,p,l; z'*\\zia);l,p,q)e'\'^\' Fccsia;^), (72) 

where z{d) is given according to ( ^OD and \\z; l,p, q) is the state (^1]). 

The u{p, q) CS representation is yet not fully specified (this could be a subject of 
separate work), except for the case oi p = 1 = q {N = 2 when it coincides with the well 
known su{l, 1) BG CS representation [§. In this case the relation ([7T|) is rewritten in the 
simpler form 

oo 

Fccs{ai,a2,'^) = Fi=o{z;<iJ) + Y,ia[ + a2)Fi{z;'^), z = aiOa- (73) 
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The BG representation is given in terms of Bargman index k, not in terms of /: 
1^) Fbg{z, k; ^f). The relation between / and k is I = ±^Ak{k - 1) + 1, or 

A; = ^(l + |/|), Z = ni-n2, (74) 
OnehasFKo(^;*) = FBG(^,fc=(l + |/|)/2;^), Fi>o(z; ^) = A;= (l + /)/2; ^) and 

Fccs{ai, Oi2\ = FBG{z,k = \l2- ^) ^Y.^(xf-^^c?t^)FBG{z, k- z = a,a2. (75) 

fc>i 

The relation (|74D stems from the definition of k by means of the Casimir operator: C2 = 
Kl — ^{K+K^ + K_K+) = k{k — 1). In the two-mode sm(1, 1) representation (H) we 
have C2 = — 1/4 + L^/4 which tell us that both I and — / lead to the same values of C2, 
that is the representations realized in the subspaces with ±Z are equivalent. However for 
the transitions between canonical CS and BG representations the sign of / is significant 
and it is taken into account in ( [73D by the identification of the standard su(l, 1) notation 
\n + k, k) of the eigenstates of once with the two-mode Fock state \n + |Z|, n) (the first 
term in the sum in (|75|) and second with \n,n+ \l\) (the second term in the sum in ([75|). 
Keeping in mind the latter identification rule we can express the two-mode canonical CS 
02) in terms of BG CS \z] k), 

\ai,a2) = \z;k= 1/2) + J2icil'''^ + af-^)\z;k), z = aia2. (76) 

k>l 

7.2 A. 2. On the overcompleteness of eigenstates of ^4^" 

Let {|^} be an overcomplete family of eigenstates of a iV non-Hermitian operator Ai, 

A.i\z) = Zi\z), 

1= j d^i{z,r)\z){z\, (77) 

where d^{z, z*) = F{\zi\, . . . , \zn\)(Pz. We note the requirement the weight function F 
not to depends on the phases of Zi. This holds for the N modes canonical CS, sp{N,R) 
and u{p, q) CS. Consider the sequence of families y^i, . . . , v^n), 

\Z ; V?!, . . . , V9„) = COS(y9„|z ] ifi, . . . , Lpn-l) + I SlU Lfn] - Z ] Lfi, . . . , ifn-l) , (78) 

where (fk,k = l,...,n, are angle parameters, n is any positive integer and z"^" is the N 
component column {zf" , . . . , z^) of eigenvalues of powers Aj" of operators Ai, 

^ 2"' I r — *9"' \ I 9^1 I — *9^ \ /'-7r\\ 

A- \ ±z = \ ±z ;y?i,...,y?„). (79) 

(fi, . . . , ipn) are superpositions of 2" states \z) with on circles of radius \zi\. If all 
(fk are integer multiple of n/2 then one gets the states \ ± z). Independent parameters 
are z,Lpi, . . . , ^pn, therefore one could also use the notation \z; Lp^"^^) (as in section 3). 

Theorem. If in equation ([T^) dn{z,z*) = -F(|-2i|, . . . , |-2Ar|) rf^f then 
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dij{z,z*)\z ■,ipi,...,ipn){(pi,...,(pn;z |, n = 0, 1, 2, . . . . (80) 

Proof. The theorem is vahd for n = hj construction. It is not difficuh to check 
directly, that it is vahd for several n > 0. Suppose now that it is valid for n — 1. Then 



that — z^" = izf ^ we obtain for the projectors in (BUI) the expression. 



we shall prove that it is valid for n as well. Indeed, using the definition ([78|) and noting 



\Z ;(fi,...,(pr,}{ipi,...,<fn;Z \=COS(fn\z ; V?! V^n-l ) ((/^i ...,(/?„„!; Z | 

+ SinV9„|-Z ■,Lpi,...,Lpn-l}{<^l...,^n-i;-Z \ + 

icos0„sinv9„|z^" V^i, . . . , v?„_i)(V9i, . . . , V9„_i; -i*^" |- 

«COS0„SinV3„| - Z ]Lpi,...,(fn~l}{Vl,---,'fn-UZ |.(81) 

We substitute this expression into equation ( [77D and then in the second and in the last 
integral change the integration variables Zi to 2;^- exp[z7r/2"~^] (rotation on angle 7r/2"~^). 
Then we note that under such rotation the eigenvalues of A^" change the sign, i.e. 

Iz"^" ^; ipi, . . . , (fn-i) I — ^ fi, • • • , V^n-i)- This yields the cancelation of the last 
two integrals and the coincidence of the first two ones in view of the rotational invariance 
of the resolution unity measure dfi{z, z*) = F{\zi\, . . . , \zN\)d'^z. We obtain 

dfi{z, Z*)\Z^" ; (fi, (Pn){(Pl, ...,(Pn;Z^"\ = 

diJ,{z, z*)\z'^" ';V9i,...,(y9„_i)(V9i,...,(y9„_i;i'^" '1 = 1- (82) 



End the proof. 

7.3 A. 3. On the uniqueness of the resolution unity measures 

dfi{zj^p,q) for u{p^q) CS 

The resolution unity measure for a given continuous family of states is generally not 
unique. It could be unique if a certain constrain is imposed on the class of admissible 
measures. For example, the requirement of invariance of the measure on the group man- 
ifold under the group action determines it uniquely |TT|. As a result the resolution unity 



measure for the group related CS is unique if it is invariant under the group action. For 
canonical SS families of noninvariant resolution unity measures have been constructed in 



ref. [HHJ. Canonical SS minimize Schrodinger uncertainty relation and can be regarded as 



group related CS for the semidirect product of 5'f/(l, 1) and Heisenberg-Weyl group |^ |. 

In this appendix we establish that the resolution unity measure for the u{p, q) CS ( pD 
is uniquely determined by the requirement the weight function F{zi, . . . ,zn^i) to be a 
smooth function of \zi\ and independent of arg^j, 'i = l,2,...,A^ — 1. Such is our weight 
function in equation (^). 

Suppose that there exists another function F'(|zi|, . . . , I-Zat-iI; q) such that the new 
measure d^' = F'd^z resolves the unity 1/ as in equation (^). Then we should have 



d^z 



F{\zp\, \zq\;l,p,q) - F'{\zi\,.. . , \zN-i\;l,p,q) \\z;l,p,q){q,p,l; z\\. (83) 
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Substituting the expansion pT|) of \\z;l,p,q) and integrating with respect to angles ipi 
argZj we obtain that the difference function 

$(ri, ra, . . . , r^_i) = F(fp, f,; l,p,q) - F'{\zi\, . . . , \zn-i \;l,P,q), 



where Vp = \zp\ = ^rf + ... + and fg = \zg\ = y^^p+i + ••• + '"at-i, should be orthogonal 
to the monomials 

r^+^..4"fi-^+\ r, = \z,\, ^ = l,...,iV-l, n, = l,2,.... 

Changing the integration variables and redenoting rf again as rj one can write this or- 
thogonality in the form 

dn... dr^-Mn, r^-iK' . . . r^'L^' = 0, (84) 



where rij = 1, 2, . . . , i = 1, . . . , — 1. Equation (|84 ) implies that $(ri, . . . , r^r-i) is 
decreasing exponentially as the total radius rf + . . . + r'^_i tends to oo. This means 
that the integral /q°° $^(iri . . . dr^^i is finite. It follows also from equation (|8^ ) that $ is 



orthogonal to any function /(ri . . .rAr_i) which admits power expansion in terms of r^, 
i = 1,2,. ..A^- 1, 

dn... drN-i^in, r7v_i)/(ri . . . rAr_i) = 0. (85) 







This implies that ^ = F — F' = almost everywhere. Indeed, if $ 7^ it must be 
nonpositive definite (in order to obey (^If)) and if $ is well behaved (it is sufficient to be 
continuous) we could find / which is negative in the domains where $ < 0. But then we 
could not maintain (^), unless F = F' almost everywhere. We suppose in (|85|) that the 
integral of the power series of / is a sum of terms of the type of (|8^) . This is ensured if $ 
is a smooth function (i.e., all derrivatives finite) of ri, . . . ,r7v-i (Our F, equation (^31), is 
such a function). In this case we can take in (^) / = $. Then we obtain that F and F' 
should coincide pointwise. Thus the resolution unity measure (^31) is unique within the 
set of smooth functions of |zi|, . . . , |2;Ar_i|. 

7.4 A. 4. Proof of the representation (|4^) of Bessel function 

In the case of g = 1 (p = — 1) and — / > our measure function F, equation (^), 
depends on ri, . . . ,rp through \z\ = + . . . + l-ZpP]^''^ = rp and it is a smooth and 

positive function of ri, ... ,rp. The measure function of ref. F' ~ fp '^P+^_ft'„/_p+i(2fp) 



is also smooth and positive [Q. Therefore the difference $ of these two functions is 
smooth and in view of (Q) and the result of the preceeding subsection they have to 
coincide pointwise. This proves formula (|46|) for Im^; = 0, Re^; > and z/ = — /— p+l = 
0,±1,.... 

Let us consider the right hand side of (^6]) as a definition of a new function F{z; u), 
z complex, u real. The integral is convergent for Re^; > and the function F{z] v) 
is evidently analytic with respect to z and v. The Bessel function K,j{z) is analytic 
and regular everywhere except of the negative half of the real line in z-plain [^. We 



proved in the above that the two analytic functions F{z; v and Ky{2z) {y = 0, ±1, . . .) 
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coincide on the positive part of the real hne of z. Then they coincide in the whole 
domain of analyticity in 2;-plain. Numerical computations show that formula (^6|) holds 
for complex v as well. In conclusion let us note that the integral in the right hand side 
of equation (^) correctly defines (under replacements \zp\ — >■ Z2) analytic 

functions F{zi, Z2; l,p, q) of two variables zi and ^2, ^^^12 > 0. At ^2 = 0, g = 1 we have 
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FIGURE CAPTIONS 



Fig. 1. Amplitude quadrature squeezing in the sp{N, R) BG CS |a; ip) and the 
superpositions |a, 0, ■0), Eq.(61), for = \ai\ — 0.05. Aqi = Aqi{f,ri,di,(j),ip), 

f = \d\, Qi = Qi ov Pi. fi = A'^qi{ri,ri,-^,(f), /s = A^pi{ri,ri, ^,0,ip), /a = 
A'^qi{ri,ri,j,0,ip), = A'^pi{4:ri,ri,j,0,tp). \d,(p) are weakly nonclassical 
states for every mode, \d, (p, ip) are strongly nonclassical. 

Fig. 2. Squared amplitude quadrature squeezing in superpositions states 
\d,(l),ip) for Tj = 0.8. AXi = AXi{r,ri,9i,(f),ip), Xt = Xi or Fj, f = 
^1 = A2X,(r„r„f,O,0) = A2F,(r„ r„ -f , 0, = A2X,(1, r„ f , 0, 0), 

53 = A2X,(1.2,r„f,0,^), 54 = 2A2p,(r„r„f,O,0) = 2A2g,(r,, n, -f , 0, t/-). 
Joint X andp (or Y and g) squeezing occurs in the interval —0.72 < t/j < —0.1. 

Fig. 3. Probabilities Pn{r,(f),ilj) to find n photons in the multimode super- 
position states \d,(f),i/j) as functions of f = \d\ for different values of and 
0- Po = Po(r, f,7r), pi = pi(f,7r, -f), = P2(r,0,7r), ps = P3(r,7r, f), 
P4 = P4(r, f , f ), P5 = P5(^, TT, -f ). 

Fig. 4. Oscillating photon number distributions p„(f, 0, -!/;) in strongly nonclassical 
states \a,(j),ip) for different values of f, and V- Pnl = Pn(0.8, 0, 7.3) (Q > 0, 
Ap>0.38, AX > 0.73), p„2 = p„(2.2,7r, -f ) (Q < 0), p„3 = p„(0.55, 5.0914, 0) 
(Q = -0, Ap>0.38, AX>1). 

Fig. 5. Nonoscillating photon number distributions in strongly nonclassical 
states |a, 0, -0) for different values of f = and 0. p^l = Pn(0.55, 2.246, 0) 
(Q < 0, Ag > 0.5, AX > 1), p„2 = p„(0.55, 2.33, 0) (Q < 0, Ag > 0.43, 
AX > 1), pnS = p„(0.55, 2.234384,0) {Q = +0, Ag > 0.5, AX > 1), p„4 = 
Poisson distribution with (a^a) = 0.685 as in p„3. 
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Fig. 1. Amplitude quadrature squeezing in 
the sp{N, R) BG CS \d; cp) and the superpo- 
sitions \d, (j), tjj), Eq.(61), for ri = \ai\ =0.05. 
Aq^ = Aq^{f,r^,(f),il}), P = \d\, qt = qt or p^. 
/i=A2g,(r„r„-f,v3), /2= A2p,(r„r„f ,0,t/'), 
/3 = A2g,(r„ r„ f , 0, V'), /4 = A2p,(4r„ r„ f , 0, V'), 
\d, Lp) are weakly nonclassical states for eve- 
ry mode, I a, ip) are strongly nonclassical. 



Fig. 2. Squared amplitude quadrature squeezing 
in superpositions states |q;, -;/') for = 0.8. 
AX, = AX,(f, r„ V'), A2 = X, or i;-, f = 

51 = A2X,(r„ r„ f , 0, V') = A2i;-(r„ r„ - f , 0, V'), 

52 = A2X,(1, r„ f , 0, V'), 53 = A2X,(1.2, r„ f , 0, V'), 
54 = 2A2p,(r„ r„ f , 0, i,)=2A\,{r,, r„ - f , 0, V')- 
Joint X and p (or y and g) squeezing occurs in 
the interval 6.4 < tp < 7.4. 
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Fig. 3. Probabilities pn{f^(j)^ip) to find n 
photons in the multimode superposition s- 
tates 1(5, if}) as functions of f = \a\ for d- 
ifferent values of (j) and ip. po = po{r,^,7r), 

Pl = pi(f,7r,-f), P2 = P2{r,0,7T), 

P3 = P3{r,7T,^), Pa = PA{r,j,j), Ps = 

P5(f, TT, -f). 



Fig. 4. Oscillating photon number distri- 
butions p)n{f'^ i'^ V') in strongly nonclassical 
states \a^(j)^il}) for different values of f, 
and ^l}. Pnl = Pn(0.8, 0, 7.3) {Q > 0, Ap > 
0.38, AX > 0.73), pn2 = p„(2.2, tt, -f ) 
(Q < 0), p„3 = p„ (0.55, 5.0914,0) {Q = 
-0, Ap > 0.38, AX > 1). 
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Fig. 5. Nonoscillating photon number 
distributions in strongly nonclassical states 
\d,(l),tjj) for different values of f = \d\, 
(f) and V'- Pnl = p„(0. 55, 2.246,0) 
(Q < 0, Ag > 0.5, AX > f), p„2 = 
p„(0.55, 2.33,0) (Q < 0, Ag > 0.43, AX > 
1), Pn3 = p„(0.55, 2.234384,0) {Q = +0, 
Ag > 0.5, AX > f), p„4 = Poisson distri- 
bution with [a^a) = 0.685 as in p„3. 
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